* JD 
gen predicted_$model_$treat=(roll1==roll2)/36
replace predicted_$model_$treat=2/36 if roll1>roll2


/*
tempvar t_pred_JD 
gen `t_pred_JD' = 1/36 
replace `t_pred_JD'=0 if roll1<roll2 
replace `t_pred_JD'=2/36 if roll1>roll2  

constraint 1 `t_pred_JD'=1 
cnsreg reports `t_pred_JD', noconstant constraint(1)  

scalar est_rmse =  e(rmse)

estat ic, df(0) n(36) 
matrix ic_results = r(S) 
scalar aic_value = ic_results[1,5]

* Store in results frame
frame change results
set obs `=_N+1'
replace treatment = $treat if _n==_N
replace model = "$model_name" if _n==_N
replace aic = aic_value  if _n==_N
replace rmse =  est_rmse if _n==_N
frame change default

*/
